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Abstract 

We study the weak convergence (in the high-frequency limit) of the 

parameter estimators of power spectrum coefficients associated with Gaus- 
sian, spherical and isotropic random fields. In particular, we introduce a 
Whittle-type approximate maximum likelihood estimator and we investi- 
gate its asympotic weak consistency and Gaussianity, in both parametric 
and semiparametric cases. 
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1 Introduction 

The purpose of this paper is to investigate the asymptotic behavior of a Whittle- 
like approximate maximum likelihood procedure for the estimation of the spec- 
tral parameters (e.g., the spectral index) of isotropic Gaussian random fields 
defined on the unit sphere In our approach we consider the expansion of 
the field into spherical harmonics, i.e. we implement a form of Fourier analysis 
on the sphere, and we implement approximate maximum likelihood estimates 
under both parametric and semiparametric assumptions on the behavior of the 
angular power spectrum. We stress that the asymptotic framework we are con- 
sidering here is rather different from usual - in particular, we assume we are 
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observing a single realization of an isotropic field, the asymptotics being with 
respect to higher and higher resolution data becoming available (i.e., higher and 
higher frequency components being observed). In some sense, then the issues 
we are considering are related to the growing area of fixed-domain asymptotics 
(see for instance [T], [13]). From the point of view of the proofs, on the other 
hand, our arguments are in some cases reminiscent of those entertained for in- 
stance by [34j . where semiparametric estimates of the long memory parameter 
for covariance stationary processes are analyzed. 

In our assumptions, we do not impose a priori a parametric model on the 
dependence structure of the random field we are analyzing; we rather impose 
various forms of regularly varying conditions, which only constrain the high- 
frequency behavior of the angular power spectrum. We are able to show consis- 
tency under the least restrictive assumptions; a Central Limit Theorem holds 
under more restrictive conditions, while asymptotic Gaussianity can be estab- 
lished under general conditions for a slightly-modified (narrow-band) procedure, 
entailing a loss of a logarithmic factor in the rate of convergence. Our analysis 
is strongly motivated by applications, especially in a Cosmological framework 
(see e.g. [S], |7]); in this area, huge datasets on isotropic, spherical random fields 
(usually assumed to be Gaussian) are currently being collected and made pub- 
licly available by celebrated satellite missions such as WMAP or Planck (see 
for instance |http://map.gsfc.nasa.gov/| ; parameter estimation of the spectral 
index and other spectral parameters has been considered by many authors (see 
for instance [13' for a review) , but no rigorous asymptotic result has so far been 
produced, to the best of our knowledge. We thus hope that the consistency and 
asymptotic Gaussianity properties we provide for our Whittle-like procedure 
may provide a contribution towards further developments. We refer also to [3], 
in , [H] J [m , [32] , [33] , [IS] for further theoretical and applied results on angular 
power spectrum estimation, in a purely nonparametric setting, and to |16) . |18j . 
[12], [H], [9], [H], [H], [19], [26] for further resufis on statistical inference for 
spherical random fields. 

The plan of the paper is as follows: in Section [2] we will recall briefly some 
well-known background material on harmonic analysis for spherical isotropic 
random fields; in Section[3]we introduce Whittle-like maximum pseudo-likelihood 
estimators for angular power spectrum coefficients based on spherical harmon- 
ics; Section |4| is devoted to consistency results, while asymptotic Gaussianity 
is considered in Section [5] In Section [6] we investigate narrow-band estimates, 
while Section [7] provides some numerical evidence to validate the findings of 
the paper. Directions for future research are discussed in Section [8| while some 
auxiliary technical results are collected in the Appendix. 
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2 Spherical Random Fields and Angular Power 
Spectrum 



In this Section, we will present some well-known background results concerning 
harmonic analysis on the sphere. We shall focus on zero-mean, isotropic Gaus- 
sian random fields T : x f2 — ^ M. It is well-known that such fields can be 
given a spectral representation such that 

Tix) = ^ ^ aimYim (x) = ^Ti (x) , (1) 

;>0 m=-l 1>Q 

aim = T (x) Ylm {x) dx . (2) 

where the set of homogenous polynomials {Yjm ■ l>Q,m — — Z, 1} represents 
an orthonormal basis for the space (S^ , dx) , the class of functions defined on 
the unitary sphere which are square-integrable with respect to the measure dx 
(see for instance [35], [H], [26], for more details, and [22], [24] for extensions). 
Note that this equality holds in both (S^ xil^dx® P) and LP' (P) senses for 
every fixed cc G S^. We recall also that for every g e SO (3) and a; G S^, a field 
T (•) is isotropic if and only if 

T{x)=T [gx) , 

where the equality holds in the sense of processes. 

An explicit form for spherical harmonics is given in spherical coordinates 
[0,7r], [0,2^) by: 



Yimi^.v) = xl ^ ^ i! , "'il ^im (cos^) e^"^^ for m > 
y 47r (t + my. 

Yim{'d,ip) = (-l)"F,,_„(t9,(^),form<0, 

Pim (cos "d) denoting the associated Legendre function; for m = 0, we have 
P/o (cosi?) = Pi (cost?) , the standard set of Legendre polynomials (see again 
[55] . [2S])- The following orthonormality property holds: 

Ylm {x) Ylm' {x) dx = 5\5m ■ 

For an isotropic Gaussian field, the spherical harmonics coefficients aim are 
Gaussian complex random variables such that 

E(a,„) = 0; 

where of course the angular power spectrum Ci fully characterizes the depen- 
dence structure under Gaussianity. Further characterizations of the spherical 
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harmonics coefficients are provided for instance by [5] , [55] ; here we simply recall 
that 

^ ^ X? , for m = , ^ for m = ±1, ±2, ±1 , 

where all these random variables are independent. Given a realization of the 
random field, an estimator of the angular power spectrum can be defined as: 

771 — — I 

the so-called empirical angular power spectrum. It is immediately seen that 

m——l \ / 

It can be convenient to write Ci = Ci/Ci . 

We shall now focus on some semiparametric models on the angular power 
spectrum; here, by semiparametric we mean that we shall assume a parametric 
form on the asymptotic behavior of C'l, but we shall refrain from a full charac- 
terization over all multipoles I. More precisely, we formulate the following: 

Condition 1 The random field T{x) is Gaussian and isotropic with angular 
power spectrum such that: 

G(i9) = G(G,ao)=G(0r"° >0, (4) 

where ao G (2, +oo) = A and for all I = 1,2, ... 

0<Ci<G{l) <C2<+00 . 

Condition A seems very mild, as it is basically requiring only some form of 
regular variation on the tail behavior of the angular power spectrum C/. For 
instance, in the CMB framework the so-called Sachs- Wolfe power spectrum (i.e. 
the leading model for fluctuations of the primordial gravitational potential) takes 
the form Q, the spectral index ao capturing the scale invariance properties of 
the field itself (ao is expected to be close to 2 from theoretical considerations, 
a prediction so far in good agreement with observations, see for instance [5] 
and [10]). For our asymptotic results below, we shall need to strengthen it 
somewhat; as we shall see. Condition [2] will turn out to be sufficient to establish 
a rate of convergence for our estimator, under Condition [3] we will be able to 
provide a Law of Large Numbers, while under Condition Hj our estimates will 
be shown to be asymptotically Gaussian and centered, thus making statistical 
inference feasible. On the other hand, in Section [6] we shall be able to provide 
narrow-band estimates with asymptotically centred limiting Gaussian law under 
Condition [2] to the price of a logarithmic term in the rate of convergence. Of 
course, the conditions below are nested, i.e. |4] implies |3| which trivially imphes 

m 
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Condition 2 Condition [7] holds and moreover, G{1) satisfies the smoothness 
condition 

G(0 = Go|l + 0(y) 

Condition 3 Condition^ holds and moreover, Gil) satisfies 

G(0 = Go|l + y+o(j) 

Condition 4 Condition^ holds with n — 0, i.e. G{1) satisfies the smoothness 
condition 

,1- 



G(0 = Go|l + o(y) 

Example 5 As straightforward examples that satisfy the previous assumptions, 
consider the rational function 

where Hi (I) and 112 (Z) are positive valued polynomials of order fc G N and 
satisfying: 

Hi (I) 

We have that 



Hi (0 ^ po+pi^ + ...pfc;'= 

112(0 Qo + Qil + ■ ■ ■ Qkl'' 



{PO 1 I I Pfc-2 1 I 111 



90 ^ 






_ gfc-2 ' 




9fc J 











Ik-l ] 1 



n, \ 90 1 J- ^ gt-a 1 I gfc-i 1 I 1 



whence 



G = ^r-(i+f^-*^U + o(l)^ . (6) 



It follows that ^ satisfies Condition ( and hence for 

^ Pk , Pk-i Qk^i 
Go = — and k — ; 

qk Pk qk 

Condition^ is satisfied when pk-i = qk-i = 0, as for instance 

1 + p ^ i-°"> r 1 



Gi^l 



-ag 



1 + 2/2 J 2 1 1 + 2/2 



1 + 1 + 2P 1 1 



or, more generally, for — ^ '^■s /o?" instance 



G, = 



1 + 2/ + 4/2 2 1 + 2^ + 4/ 
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3 A Whittle-like approximation to the likelihood 
function 



Our aim in tliis Section is to discuss heuristically a Whittle-like approximation 
for the log-likelihood of isotropic spherical Gaussian fields, and to derive the cor- 
responding estimator. Assume that the triangular array m = 
I = 1,2, ...L is evaluated from the observed field {T{x)} , by means of ^. Our 
motivating rationale is the idea that a set of harmonic components up to mul- 
tipole L can be reconstructed without observational noise or numerical error, 
whereas the following are simply discarded; this is clearly a simplified picture, 
but we believe it provides an accurate approximation to many current experi- 
mental set-ups. Of course, L grows larger when more sophisticated experiments 
are run (L can be considered in the order of 500/600 for data collected from 
WMAP and 1500/2000 for those from Planck). It is readily seen from ^ that 
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where the variables {aio, -\/2 Re {an} , a/2 Im {an} , Re {an} , \/2Im{an}} 
are i.i.d. Gaussian variables with law JV{0, Ci), see [2]. The likelihood function 
can then be written down as 

-2 \ogCi{d; {a,™}L=-i) = (2/ + 1) log(27r) + [21 + 1) log Cm + {21 + 1)^ , 



or equivalently 



const + {21 + 1)7^ - {21 + 1) log 



'GW ' ' °Ci{d) 

Clearly this landscape is overly simplified, for instance, due to numerical 
errors and aliasing effects the expected value E |a/m|^ may not be exactly equal 
to the population model Ci {•&) ; however in Conditions [T] and following we are 
allowing the two to differ to various degrees, and we expect this to cover to some 
of effect these experimental features that we are neglecting. More generally, 
rather than a sharp cutoff at L, a smooth transition towards noisier frequencies 
would represent more efficiently actual experimental circumstances, but we view 
our results here as a first step amenable of further developments. 

An alternative heuristics for our framework can be introduced considering 
that for I = 1,2,...,L, the following Fourier components can be observed on a 
discrete grid of points {xi, xk} 

~fi = {Ti{xi), ...,Ti{xk), ...,Ti{xk)} ■ 

To simplify our discussion, we shall also pretend that {xi, ...,xk} form a set 
of approximate cuhature points with constant cuhature weights Xk = At: / K (see 
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for instance [5S], [12]), so that we have 



X] -^^'"11 {xk}Yijn2 (xk) ^ , for Z = 1, 2, L . 



K 



Again this landscape is overly simplified, for instance, cubature weights on the 
sphere are known not to be constant, but their variation is usually considered 
numerically negligible. The frequency components T; are well-known to be in- 
dependent and we can hence write down the likelihood function as 



1^ 



1=1 



where 



21 + 1 
An 



CiPi[{xj,Xk)) 



The matrix 51/ can be (approximately) decomposed as follows: 



Yi,^i{xi) Yi,^i+i{xi) 
yi.-i{x2) 

Yi-i{xK-l) 

Yi,^i{xk) Yi,^i+i{xk) 



Yu^i{xi) Yi^iixi) 
Y1AX2) 

Yi^l{xK-l) 

Yij^iixK) Yi^i{xk) 



X— C;/2/+lX- 

47r 



In fact 



Hence 



Yi,^i{xi) y,,-,(x2) 

Yi,^i^i{xi) 

F^;-i(a;i) _ ... 
Yi,i{x^) Yu{x2) 

=:yi X Ci{d)hi+ixy; 



Yi^^i{xK-i) Yi.-i{xk) 
Yi,^i+i{xk) 



Yi^l{xK-l) 



ytyi ^ hi+i and det{f^,} ~ Cf+'i^) 



Yjj_i{xk) 

Yi,i{xk) 



-2\og Cii^-^i) ^K+{2l + l)\ogCi{d) + [f[yi X Ci\d)hi+i X y*i^i] 
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Now 



y0i = 



Y 



1,1-1 



Yi-i{xK-i) Yi^_i{xk) 
Yi-i+i{xk) 



Y 



Yu{xi) Yu{x2) ... Yi^i(xK^ 

O-hnYlmixi) 
Dm aimYlrn{x2) 



0-lmYlm{xK-l) 
Z^m aimYlmixK) 



J2,n, O'lrni J2kYl-niliXk)Y l-l{xk) 
Emi "'"il J2k^l™'li^k)Yl-l+i{xk) 

Emi a;mi J2k^lnii{Xk)Yj^l-i{Xk) 
Emi aim, T,k'^lrnAxk)Yl^l{xk) 



1,1-1 



Yu{xk) 





ai,-i 




ai,-i+i 


~ V 47r 






a-Li-i 




ai,i 



whence 



A-K 1 



K am 

As before, we can then conclude heuristically that 



2\og Ci{^;l^i) ~ const + (2/ + 1) 



Ci 



cm 



{21 + 1) log 



Ci 



(7) 



Again we stress that for a general spherical random field with an infinite-terms 
expansion such as ([I]) the relationship ([t]) cannot hold exactly; indeed, precise 
cubature formulae can be established only for finite order spherical harmon- 
ics. In general, this may introduce some numerical error: as mentioned before, 
however, we pretend in this paper that such correction factors are covered by 
Conditions [T]|4j In other words, we envisage a situation where data analysis is 
carried over on multipoles / where numerical errors are of smaller order and the 
approximation Q holds for the expected variance of the sample {aim}. 



4 Asymptotic results: Consistency 

As motivated in the Introduction, in this paper we shall not assume we have 
actually available a fully parametric model for the angular power spectrum. 
Instead, the idea will be to use an approximate maximum likelihood estimator, 
which shall exploit the asymptotic approximation provided by Condition [ij i.e. 
Ci ~ Gl^". In view of the discussion in the previous Section, the following 
Definition seems rather natural: 
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Definition 6 The Spherical Whittle estimator for the parameters (ao,Go) is 
provided by 



(S.,G.):=arg^mm^g|(2/ + l)^-(2/ + l)log^| . (8) 

Remark 7 For general parametric models Ci — Ci{'&), the Spherical Whittle 
estimator for a parameter g O C W can be obviously defined as 

Remark 8 To ensure that the estimator exists, as usual we shall assume through- 
out this paper that the parameter space is a compact subset of AxT gM.'^; more 
precisely we take a £ A — [01,02] , 2 < ai < 02 < 00, and G £T — [71,72] , 
< 7i < 72 < 00. This is little more than a formal requirement that is standard 
in the literature on (pseudo-)maximum likelihood estimation. 

We can rewrite in a more transparent form the previous estimator follow- 
ing an argument analogous to [34], i.e. "concentrating out" the parameter G. 
Indeed, the previous minimization problem is equivalent to let us consider 

{aL,GL) : = argmin7?.L(G', a) 

a.G 

TZL{G,a) : = ^(2^ + 1)^ + ^(2/ + 1) log G + ^(2Z + 1) log T". 
i=i ;=i ;=i 

It is readily seen that: 

1=1 1=1 



so that 



57^L(G,a) , _ ^ , 1 , Ci 



dG 







G = G{a) := V (2/ + 
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Moreover, we have 
9G2 



1=1 



1 

G2 



G2 

L 



9G2 



1=1 

L 



G=G{a) 



G(a) Z^(=il^' + J-ii^ i=i 



1 

^^(2? + l)>0. 

G{a)\ 1=1 



Considering, moreover, that on the boundary of A x F, TZl{G, a) — +oo, G — )■ 
TZL{G,a) has a unique mininimn over its domain on G = G(a). Thus we can 
consider: 



ul = arg min _Ri(a), 

a 

nLiG{a),a) 



(10) 



Er=i(2? + 1) 
= (logG(a) ^ y(2/ + l)logn. 

V Ef=i(2^+i)tr ; 

Our first resuh is the foUowing 

Theorem 9 Under Condition^ as L oo we have (compare ^) 
moreover, under Condition^ 

Gl -^p Go . 

Proof. To estabhsh consistency, we shah resort to a technique developed 
by [H] and [31]. In particular, let us now write 



G{a) G(ao) (a - uq) 

log TTTT - log 



G(a) 



where 



log 



G(a) 

G(a) 
G{ao) 



Z^;=i(.^' + Ij (^-^ ' 

logl^ y(2/ + i)r-"4 . 
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so that 



Ai?L(a,ao) = ULia,ao) - TL{a, ao), 



The proof is then completed with the aid of the auxihary Lemmas [TT] [12] 
that we shaU discuss below. Indeed 

Pr{|aL - aol > e} < Pr i inf Ai?L(a,ao)<0 

|^|a-Q!o|>e 



< Pr inf [UL{a,ao)-TLia,ao)]<0 

l\a-ao\>6 

For ao — a < 2 the previous probability is bounded by, for any S > 



and 



<Pr<{ inf UL{a,aa) <S} +Pr{ sup TLia,ao)>0} , 

Jq-qo|>£ J l|Q-Qol>e I 



lim Pr <^ sup Tl (a, ao) > ^ = 



|a-ao|>e 



from Lemma 12 while from Lemma 11 there exist 6^ ^ {1 + e/2) — log(l + 



e/2) - 1 > such that 

lim Pr j inf Ul (a, ao)<sA ^0 . 

i— i-oo — Qo|>e J 

For ao — a = 2 or ao — a > 2 the same result is obtained by dividing Ai?L(a, ao) 
by, respectively log log L or log L and then resorting again to Lemmas [Tl} [T2| 
Now note that 





-1) 


Go 






'1) 


1 




Ef=i(2H 


-1) 


1 






f 1) 



Er=i(2^ + i) 

Ci ' 



X^(2Z + i)Gor("°-"-) I - I = 



5^(2/ + l)Gor("°-3.) J / - 1 ] + (i - /(-o-s.)^ 



Go/-"« 
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Clearly: 

G(Si) - Go 



< 



1 

Ef=i(2/ + i)tr' 



Ci 



Go 

Ef=i(2/ + i)fet 

\Ga\ + \Gb\, 



Gol-'-o 

Y^{21 + 1) (l-^o'o-^^) _ 



so that 



Pr(|G(aL)-Go| >£) < Pr (|Ga| > |) +Pr(|GB| > 
Observe that: 



Pr{|G.|>|}_ 
+ Pri \ao-aL\ > 



:r !■ < Pr 

1 



\GA > ^ 



n 



|Q!o - ail < 2 



< Pr ■ 



1 

< 



1 



Er=i(2^+i)tr 
1 



^(2i + i)Go;^/' 



Ci 



GoZ-«o 



Go/-"" 



< 



G 



1 



^E"=i(2/ + l)S 

G L"/6 



^(2^ + l)Gol^/^E 
-1 

X^(2; + i)Go|i/3ri/^ + oz.(i) 
1 

ol{1) = ol{1) . 



> £ 



ol{1) 



ol{1) 



^E,=i(2« + 1) 

As far as the second term is concerned, we have, for a suitably small i5 > 0: 



Pr 



{\G, 



> 



Pr 



Pr 



{[\G, 
{[\Ge 



> 



> 



n [log I {ao - Sl)] < 5) + Pr (log I (ccq - Si) 

r)[\ogl{ao-aL)] <S^ +ol{1) ■ 



and using |e ^ — 1| < a; for < a; < 1, we obtain 



/-(ao-Si,) _ ^ 



|exp(-log/ (ao - Sl) - 1| < log? |ao - Sl| , 



Pr {[\Gb\ > |] n [log/(ao - a^)] < 5) 
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< Pr 



< 



< 



> - n [logZ(ao - oil)] < S 



Ef=i(2; + i)tt ^ 



where we have used 



E|ao - SlI < {e |ao - = O(^) 



which under Condition [2] wiU be estabUshed in the proof of Theorem 16 



The first auxiliary result we shall need concerns G, G and their fc-th order 
derivatives Gk,Gk, i-e. 



Gk («) 



Gk (a) 



E,t,(2l + 1) S 

L 



1 



E(2Z + l)(logO'^^ ,fc-0,l,2,.. 



where Go (a) = G{a) and Go (a) = G(a) defined as above. 

Lemma 10 Under Condition^ for all 2 > ao — a > £ > 0, as L — > oo, we 
have 

Gfc(a) 



sup 



log 



Gk{a) 



On the other hand, if ao — a > 2, 



sup 



lot 



Gfe(a) 



Gfe(a) 



Proof. Let us first focus on the case where a — ag > —2 . For clarity of 
exposition, we start from a simplified parametric version of Condition [l] i.e. we 
assume that we have exactly 



= G(Go,ao) -Gor"°. 



(13) 
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Let us write first 

Gfc(a) 



Gfc(a) 



1 



Ef=i(2/ + 1) (logo' 1% - ElA^i + 1) (logo' ^ 



Eti(2^ + l)(log/) 



Eti(2/ + 1) (logo' GoP-"« {^Jb^ - l} 



Eti(2Z + l) (logO'Go^"- 



Fixed < /3 < ^, we have, for all I : 



Pr 



Ef=i(2Z + l)GoP-"« (logo' {gS^ - l} 



Eti(2/ + l)(logO'GoZ"-« 



< Pr L'' 



Ef=iy(2^Ti)(iogO'^"-"'' 



(2^ + 1) (log o';"-"« 



sup, ^{21 + 1) 



- 1 



< Pr l^sup \/(2? + l) 
because 



Ci 



^;3 Eiiiy(2m)(logO'^°-"° ^ ^ X/^+l+«-^o log^L 



Ef.i(2? + l)(log0'?"-"'' 



= Gi'^-1/2 = o(l) 



Now 



Pr I sup v/(2/ + 1) 



Gz 



< L max Pr <^ {21 + 1) 



GnZ-"o 



G, 



GoZ-"o 



- 1 



and 



Pr<^ V(2; + l) 



G, 



Go/-"" 



- 1 



E 



> ^,L'^ \ < C- 



7(21+1) 



= 0(L- 



uniformly in I, see for instance such that Af > 1//3. Hence 

G, 



Pr <^ sup ^/{2l + 1) 



GoZ-"° 



- 1 
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For the general semiparametric case, the only difference is to be found in the 
expressions for ECi, which under Condition [2| becomes 

ECi^Gol-"" (i + o(ri)), 

where the bound O {l^^) is uniform over a by Assumption. As before, we hence 
obtain 

Gkia) 



1 



Eti(2/ + 1) (logo' 1% - Ef=i(2^ + 1) (logo* 



Eti(2^ + i)(iog0 



Eti(2^ + l)(logO'GoP-"°{ 



_Ci 



Go/-°o Goi-°o ' Gol 



EGi _|_ ECi 



Eti(2/ + l)(logO'Go;"-"o 



Eti(2/ + l)(logO'GoP-""{ 



a — an } Ci EC; 



Ef=i(2/ + l)(logO'Go/"-"« 

^ Ef=i(2/ + i) {\ogif G^r-^^ {o{\)] 

Eti(2Z + l)(log/)'=Go;"-"o 

The second summand is immediately observed to be O(j^). By the same argu- 
ment as before, for < /3 < |, we have, for all I : 



Pr 



Ef=i(2/ + l)(logO'Go;"-""{ 



EG; 



Ci 



Ef=i(2Z + l) (logO^Go^"-"" 



>5, 



7(21 + 1) (log o'P 



E,(2Z + l)(logO'Z"- 



sup, ^{21 + 1) 


C, ECi 


Goi"°o Go!"°o 




LP 



EG, 



< -pv/(2m)-^^^_ 



EC, 



< PrJsupV(2; + l)|l + 0(y)| 



- 1 

Ci 
EG 



> 



The rest of the proof is analogous to the argument we provided before, and 
hence omitted. 

For the case where — a > 2, it suffices to note that 



Gfe(«) 1 



|^^-t5ITT)Er=i(2^ + l)(logO 



> , with probability 1 , 
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and 

^ Sd^i:t.(2/ + I)(log0'-^{1 + 0(})} ^ 
G^(-) ^^-feTyEti(2^ + l)(log/f ^ 

■ 

We are now in the position to establish the asymptotic behavior of Ul (a, ao) 



in (111, for which we have the following 



Lemma 11 For all 2>ao — a>e>0, we have that 



lim I- i"'""^ E(2/ + l)log/ + log-^| 



= lim 



log I V(2? + l)r-- I - V(2/ + 1) log I 



= (1 + (a - ao)/2) - log(l + {a - ao)/2) - 1 ><5, > . 
Moreover, if — a =2, 

lim -4- 1 - - V + 1) i«g ^ + I - 1 > ' 

L^^loglogL[ Ef=i(2^ + 1) G(ao) J 

and /or ao — a > 2, 

lim 1"'"°^ V(2/ + l)log^ + log^Uao-a-2>0. 

L^<x,logL| Ef=i(2^ + l)tt ^G(ao)j 

Proof. Note first that 

lim ^^^o) v-^2Z + 1) log / = (a - ao) , (14) 

^-ooiogixEf=i(2« + l)tt 



and 



hm {logLj-^logj ] V(2^ + ir-"4 =(«-ao) , (15) 



with 



log I ^-^ V(2Z + i)r-- 1 - V(2; + 1) log z > 

for all L > 0, because by Jensen's inequality 

logEp,, {a:"-""} -Ep„, log{a;"-"°} > 
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where 

P / n 2; + l 

"''^=''*^ = " = SLP^TT)- 

To prove that the inequaUty is strict also in the hmit, consider first the case 
a — ao > —2 

log J y{2l + I 

= log I J-ao^f i'(2!+^l) ^^^^ ^ 1)^"-"" |-log(l+(a-ao)/2)+(a-ao) log L, 



where 

L 



(l + (a-ao)/2) ^(2; + i);«-ao_i = ^^i) , whence 



log/ + V(2/ + i)r-"4 = 0,(1). 



Thus 



■ log(l + {a- ao)/2) + {a- ao) log L + - Jl^ ^ ^ V(2« + 1) log / + 0,(1) 



Now 



(a - ao) ^ 



Li=l(2' + 1) ; = i ^ 

1=1 

. X , L dx (a — an) 
= 2(a-ao)y^ Z^°^xL-^-^+^^(^) 

= -2(a-ao)/ xloga^da; — ^ + 0,(1) = ol(1) 

Jo ^ 
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because 



/ X log xdx 
-I 

We have hence proved that 

1 



log a; 



Jo 2 X 4 



log • 



J2i2l + l)l 



a-ao 



Eti(2/+i) 

= (1 + (a - ao) /2) - log(l + (a - ao) /2) - 1 + 0^(1) > 

for all |a — ao| > £, a — ao > —2. 

Consider now the case ao — a > 2. We can rewrite: 



^(2Z + l)log/ 



G{a) 



. y(2Z + l)log/ + log^. , 



= log • 



Eti(2? + l)t^ 



^(2Z + 1)Z° 



Eti(2^ + i)tr 



^(2/ + l)log/ 



= (ao - a) log L 
For the term A^: 

L 



log Et 1 (2Z + i)r log Ef=i (2Z + 1) 



log; 



(ao-a)logL (ao-a)logX ' Ef=i(2« + 1) ^""^^ ^^^og^ 

= (ao - a) logL [Al + + Cz,] . 



^(2; + l)l-(-o-a) ^ ^J2l^-(ao-a) ^ o^.-^^,-., (1) C> 1 



(=1 



1=1 



because Ef^i 1^ a convergent series when the exponent 1 — (ao — a) < 

-1; for 1 - (ao-a) = -1, we have |Ef=i ^^"^"°""Vlogi>} 1 and the 
argument is analogous. Therefore 



(ao - a)logL x [Al] = 
As far as Bl is concerned, we have: 

L 



O (log log L) for ao — a = 2 
0(1) for ao - a > 2 



log^(2Z + l) = log L^ + o (log L) = 
1=1 

= 21ogi + o(logL), 



so that: 



lim Bl 

L— >oo 



(ao - a) ' 
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finally, simple manipulations and standard properties of the logarithm (which 
is a slowly varying function, compare yield 



lim Cl = lim 



Summing up, we obtain: 



El 



- — y:(2?+i 



logL 



= 1 



lim {(ao - a) log L[Bl +Cl]} = 



for aQ — a = 2 
{ao — a) — 2 > for ao — a > 2 



and the claimed result follows. ■ 

In [34' a related computation was given for approximate Whittle estimates 
on stationary long memory processes in dimension d = 1, i.e. the limiting lower 
bound turned out to be (1 + (a — ao)) — log(l + (a — ao)) — 1 + ol(1) > S^- 
In view of this, we conjecture that for general c?-dimensional spheres the lower 
bound will take the form 

(l+^°^)-log(l+(°^)-l + o,(l)>^.. 

Now we look at (a, ao), for which we provide the following 
Lemma 12 Under Condition^ as L oo, we have 

sup ITl (a,ao)| = Op(l) , /or ao - a < 2 , 

a 

sup iTi, (a,ao)| = Op(l) , /or ao - a > 2 . 



Proof. For ao — a < 2, consider first 



Gjao) 
G{ao) 



1 



1 



Er=i(2^ + i) 

where we have easily, as L — > oo. 



Gi 



Go/-"° 



Gjao) 
G(ao) 



1 



1 — 

■2/ + l)^ Gol-"" 



Er=i(2^+i) 



Var 



I Gjao)] 
\G{ao) j 



^ 1 
E(2/ + l)2 = 0( 



1=1 



Eti(2^ + 1) 
whence by Slutzky's lemma 

On the other hand, in view of Lemma [TO) we have that: 

G(a 



sup 



log 



G(a) 



Op(l), 



whence the result follows easily. The proof for ao — a > 2 is immediate. 
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5 Asymptotic properties: The Central Limit The- 
orem 

For the Central Limit Theorem, we shaU exploit some classical argument on 
asymptotic Gaussianity for extremum estimates, as recalled for instance by j30j . 

Theorem 13 (Newey Mc Fadden CLT) Assume that the extremum estimate 
OLM '■= argminctgyi R{a) satisfies the following conditions, as N ^ oo: 

^ r 

1. an — > ao; 

2. ao ^ dA 

3. Rn {a) e 6*2 (/^ (ao)) , for all N; 
4- i?7v (a) = a.s.; 

5. ViVi?^(a) AaA(0,S); 

6. there exists F — F {a) , F (ao) ^ such that 

sup ||i?^(a)-F(a)|l -^0 , 

Q:|a — ao|<e 

The ^ 



N{a-ao)^Af (o, (E/F (ao))') 



The proof is standard textbook material - here, we simply recall a few steps 
to fix some notation we shall use later. Conditions 2, 3, 4 in this Theorem imply 
that R'pf (a) = a.s.. We can hence use a Taylor argument to expand i?^ (a) 
around the point ao: 

= R'l^ (a) = i?^ (ao) + (S - ao) R'n (a) + o ((S - ao)') , 

so that: 



, NR'^ (ao) 

TV (a - ao) = , 

(a) 

where a is a mean value that lies on the line segment between a and ao- Because 

p 

a — > ao and Condition 4, with probability approaching 1, we have 

|i?^(a)-F(ao)| < |i?^(a)-F(a)| + |F(a)-F(ao)| 

< sup|i?^(a)-F(a)| + |F(a)-F(ao)| ^0. 

p 1 P 1 

Hence F (a) — > F (ao) and because F (ao) ^ 0, F (a) — > F (ao) and the 
result follows easily by the Continuous Mapping Theorem. 

We now introduce some auxiliary results whose proof is deferred to the Ap- 
pendix. 
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Proposition 14 Let 



Zl (s) : 
Then, for s G 



L L / L \ ' 

^ {21 + 1) + 1) (log 0' - E + 1) ''''' ^ 



Moreover, let Li ~ L-{1 — g (L)), where < g (L) < I is such that limL^ao 9 {L) - 
0. // 

L L / L \ 2 

l=Li l=Li \l=Li ) 



(17) 



have 



/™ 4fi+.l («) = ^ (5) , (18) 



where 



12 8 3 



2^ 

Remark 15 It is obvious that, if s ~0, 

Ko ^ K is)l^, = ^ . (19) 

Our purpose below is to show that the conditions provided in Theorem [T3| 
in are satisfied by and Rl in the Spherical Whittle estimates, so that the 
following result will be established. 

Theorem 16 Let oil = argmiuagyi RL{a) defined as in ([i 

a) Under Condition^we have that 

|e (Si — Q!o)^| = Q( ^ ) , whence (Sl — ao) = Op{ ^ ) , as L —>■ oo . 

(20) 

b) Under Condition^ we have that 
^ (Sl - ao) — >p -K . (21) 



41ogL 

c) Under Condition^ we have that 



4^(aL-ao) AaA(0,1) . (22) 
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Proof. We note first that under Condition |4] ([22| ) is an immediate con- 
sequence of (21 1; on the other hand, the proof of (|20|) follows on exactly the 



same lines as (21 ), the only difference here being that the asymptotic bias term 



cannot be given an analytic expression but only bounded. It is then sufficient 



to establish (21 ), as we shall do below. 



Following the notation introduced above, for each L there exist G (ao — a, ao 
such that, with probability one: 



{oil - ao) = - 



where 5*/, (a) is the score function corresponding to Ri^{a), given by: 



da Gia) 



1 ^ 
y(2Z + l)log/ , 



and 



da 



SLia) 



da 



G2{a)G{a) - Glia) 
G^a) 



ElA^i + i)^ {Ef=i(2; + 1) A} - {Etii^i + i)(iog/)^}' 

{Ef=i(2^ + i)ia}' 

where G{a), Gi{a), 6*2 (a) are respectively the estimate of G and its first and 
second derivatives, as in LemmafTOl By direct substitution we have immediately: 



SL(a) 



1 ^ 

V(2Z + l)logi 



Ci 



Now, 



SL[ao) 



Go 
Go 



G(a)l- 



Ci 



G(ao) 
Go 



G(ao) 



where 



SL{ao) 

Go 
G(ao) 



y(2/ + l)log/. 



Ef=i(2^+i) 

1 + Op(l) , as L ^ oo , 



Ci 
Go/-"° 



- 1 



and 
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in view of Lemma [Tol Also 

L 



E5i(ao) = —r^ y(2/ + l)log/(— ll 



and 



lim 2LVar{S'L(ao)} = 1 (23) 



In fact, we have: 

Var{SL{ao)]=Vi + V2 + Vs 



where 



2^(2^ + 1) (logo' 



1 ^ 2 r 

Now because 

= ; (24) 



Cow 



Q G(ao) 
Ci ' Go 



we have 



Er=i(2Z' + l),'tl 
2 

Er=i(2/ + 1) 



Far{5i(ao)} 



E(2r + l)Go.(^,^'' 



Ci ' Gr 



(25) 
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2 _ 1 
Z^T ~ 2L2 



—3 [ ^ (2/ + 1) ^ (2; + 1) (log if - (2^ + 1) log l] 
1)1 \;=i 1=1 \i=i ) 



by using (16) and (30) with s = to obtain (23). In order to establish the 
Central Limit Theorem, it is sufficient to perform a careful analysis of fourth- 
order cumulants (note our statistics belong to the second-order Wiener chaos 
with respect to a Gaussian white noise random measure). Write: 



where 



Ai = (2; + i)iogr 

Bi = (2; + l)logr 
In the Appendix, we show that 



1 



Gi 

GLjao) 
Go 



(26) 
(27) 



l^cumlj2^A,+Bi,),Y,{Ai,+Bi,),Yl (^'3 + Bi.. ) , ^ (A, + B,, )\ = Ol( 

\ ll 1-2 Is ^4 y 

whence the Central Limit theorem follows easily from results in [31]. Indeed, 
using recent results from the latter authors a stronger result follows, that is 

dTy(5]X,^L,^) = 0(7) , ^ = AA(0,1) , 



log^L 

L2 > 



1=1 



where (iTV'(VF, V) denotes the total variation distance between the random vari- 
ables W,V, i.e. 



drviW, V) = sup |Pr {W e B} - Pi {V e B}\ , any Borel set B . 



Also 



L 



-ESLiao) 



''Ef=i(2/ + i)tt 



{21 + 1) Yogi 



0(1) 



log L . rsi/ + n ^ I log L 

— > — K , as L — > cx) . 

Let us now focus on the second order derivative. From consistency, it is sufficient 
to focus on |a — ao| < 2; here we can apply again Lemma 10 replacing the 



24 



random quantities Gk{cy) with the corresponding deterministic Gk{ot) values, to 
obtain 

G2{a)G{a)-Gl[a) 
QL[a) = OHa) ' ' 

uniformly over a. It is convenient to write 

G2{a)G{a) ~ Gl{a) _ Ql^"'{a) 



G2(c 



Let us start by studying Q'j^^{a). We have, by using (301 with 5 = and 
s = a — ao- 



^den 



= 2 + 

V(i + ^) 

Consider now QJ^"''{a), where we have: 

L2("-"o) 



^2 

^4+2(Q-ao) 



^(2/ + !)— log^Z E(2^ + l) 



7-ao \ I-ao 



+ 1) — log' i ^(21 + D-p^ - + 



w=i 



^^ = l 



^2 



1 



4(1 + ^. 



by using (16), s = a — uq Combining all terms, we find that, uniformly over a 
G^ . 1„..4 +ol(1) 



1 



Gn 



Finally, from the consistency result 



4(1 + ^, 



+ Oi(l) . 



1 



Thus by Slutzky Theorem we have: 

4 QdaL) ' 



1 , QL{aL) 
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as claimed. ■ 

In the Appendix we describe in details the result concerning the analysis of 
fourth-order cumulants. 



Remark 17 In the statement of the previous Theorem, we decided to report 
normalization factors in the neatest possible form. A careful inspection of the 



proofs reveals however that the asymptotic result in (21) and KM) can he im- 

L 1^1 = 1 logL 



proved in finite samples introducing a correction factor cl ~ ^a=i ~^ 



as L oo, as follows 



41ogL X cl 

under condition Q), and 

V2L 



Ax cl 



under condition Note that c^ < 1 for all finite L, whence the asymptotic 
bias and variance are slightly underestimated in Theorem \16[ For instance, 
the correction factors for L = 1000, 2000, 4000 are, respectively, ciooo — 0.86, 
C2000 — 0.87, and C4000 — 0.88. 

The previous result provides a sharp rate of convergence for the spherical 
Whittle estimator. However in the general case the asymptotic bias term —^/2k 
is unknown, which makes inference unfeasible. To address these issues, we shall 
consider in the next Section an alternative narrow-band estimator (compare [34]) 
which achieves an unbiased limiting distribution, to the price of a log factor in 
the rate of convergence. 



6 Narrow- band estimates 

In the previous Section, we have shown that under Conditions [2j [3] it is possible 
to establish a rate of convergence for the spherical Whittle estimates; however, 
due to the presence of an asymptotic bias term, statistical inference turned 
out to be unfeasible. The purpose of this Section is to propose a narrow band 
estimator allowing for feasible inference under broad circumstances. We start 
from the following 

Definition 18 The Narrow-Band Spherical Whittle estimator for the parame- 
ters = {a,G) is provided by 

(SL;L,,Gi;Lj :=argminE |(2/ + l)^-(2Z + l)lor ^' ' 
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or equivalently 

aL-Li = arg niini?L.Lj (a, G'(q:)), (28) 

a 

Rl^lM.G) = (logG^;^,M- " V(2/ + l)log^| , 

where Li < L is chosen such that 

i — Li— >oo . — = 1 + 0(- -) , as L —> oo . 

' Li ^ogL 

For definiteness, we shall take 

Li= L{1- g (L)) , where g (L) 



logL 



pretending for notational simplicity that Li so defined is an integer (if this isn't 
the case, modified arguments taking integer parts are completely trivial). 

Theorem 19 Let Sl;Li defined as in (28). Then under Condition^we have 

^^^^(a.;.,-ao) Aaa(o,i) . 

Proof. The proof of the consistency for aL;Li can be carried out analogously 
to the argument provided in Section |4j and hence is omitted for brevity's sake. 
The proof for the Central Limit Theorem can also be carried along the same 
lines as done earlier, noting in particular that for the form ^ of C; 

ESl-mM = r ! ^ V (2^ + 1) {log/} ' 



EU(2« + l).f^; [Go/-- Go 



L 



{21 



I 



EtL,(2? + l), 



logLi (\ogLi\ logLi 

K- 



M=o(^), 



and 



= O 



Li \ Li J Li 



L-g{L)^E[SL;LAc^o)] 
log Li\ L 



^1 J log2 L 



0(f)0 i^Uop^ =0.(1) 
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On the other hand 



Var {SL-Liiao)} 



1 



Ef=L,(2^ + l) 
2 



rVarl J2i2l + l){\ogl} 



Ci Gl:Li (a) 



J=Li 



Gq 



Er=L,(2^ + i) 

2 



L L / L \ • 

^ {21 + 1) ^ {21 + 1) {log^ /} - ^ (2/ + 1) {log/} 



Ef=L,(2^ + l) 



3^L;g(L) (0) 



by using (24 1 and (25) and foUowing the notation of Proposition 14 with s — 0. 
Proposition [14] leads to: 



while 

L 



-I 3 



^(2/ + l) 



= {L' ~ L?)-^ = 8i« ( g (L) - '-g' {L) ) = SL^g' {L)+o, {L^g^ {L)) 



By substituting these results in ( 29 ) we obtain 

9{L) 



Var {S'L;Li(ao)} 



Rewrite now the term (5lil(q;) as 



QLi-L{a) 



where we have:\ 



1 



127.2 12^2 log (2,) 



Qnum / „ \ 



^2 



(Eti (2/ + 1) 



r2'L,g(L) (s) , 



EtL,(2^ + l). 



E(2' + ir 



where s = a — ao- 
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From ( 30 1 , we have 



Qf"L(«) - Gl ) ^^+o^(i) 



+ (i-(l-.g(L))f 



«l(1) 



Consider now Q2^"l{a), where we have: 
Combining the two results, we obtain: 

L— *-oo 4 

Finally, from the consistency results, we have: 

12 

Qh^-Lia) -^p 1 . 



9' (L) 

The analysis of fourth-order moments is exactly the same as in the previous 
Section, and the result follows accordingly. ■ 

Remark 20 It should be noted that an asymptotic unbiased estimator is ob- 
tained with the loss of only a logarithmic term to the power 3/2 in the rate of 
convergence. 

7 Numerical Results 

In this Section we present some numerical evidence to support the asymptotic 
results provided earlier. More precisely, using the statistical software R, for 
given fixed values of L, ao and Go and the alternative conditions discussed in 
the previous Section, we sample random values for the angular power spectra Ci 
and we implement standard and narrow-band estimates. We start by analyzing 
the simplest model, i.e. the one corresponding to Condition |4] Here we fixed 
Go = 2. In Figure 1 we report the distribution of Si, — ao normalized by a factor 
y/2L/4:. In Table 1, we report instead the sample frequencies corresponding to 
the quantiles q = 0.05, 0.25, 0.50.0.75, 0.95 for a TV (0, 1) distribution. 
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L=ZaDO 



L=3DDD 



L=4000 




alphaO = 2 



■i ; 

alphaO = 3 



-4 -i c ; 

alphaO = 4 




Figure 1: Distribution of normalized (S^ — ao) by varying L and aQ, under 
Condition IH 
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Sample 


Frequencies 








L 


-1.96 


-1 


-0.68 





0.68 


1 


1.96 






2000 


4 


19.2 


29.2 


48.6 


22.8 


14.2 


4 




2 


3000 


4.5 


18.4 


26.8 


51 


23.33 


14.36 


3.6 






4000 


4.4 


17.7 


25.2 


49.1 


23.43 


13.87 


4.8 






2000 


4.4 


19.2 


29.2 


51.5 


24.03 


15.17 


3.6 


ao 


3 


3000 


4.3 


18.4 


26.8 


48.9 


23.2 


13.43 


3.8 






4000 


4.2 


17.9 


26.4 


50.8 


23.07 


14.13 


3.7 






2000 


4.4 


21.6 


30.2 


50.9 


22.73 


14.94 


5.5 




4 


3000 


4.2 


21.2 


29.8 


50.4 


25.07 


15.87 


4.3 






4000 


4.2 


17.9 


27.1 


50.4 


22.7 


13.73 


4.2 



Table 1: Quantiles of L/4\/21ogL (S^ — ao) 



Table 2 provides the results for the classical Shapiro- Wilk Gaussianity Test 
performed on simulations obtained by varying ao and the number of multipoles 
L. Asymptotic Gaussianity is clearly supported. 







Shapiro 


-Wilk Test 




L 


W 


p — value 






2000 


0.9976 


0.685 




2 


3000 


0.9978 


0.667 






4000 


0.9983 


0.373 






2000 


0.9976 


0.691 


ao 


3 


3000 


0.9980 


0.842 






4000 


0.9985 


0.945 






2000 


0.9987 


0.670 




4 


3000 


0.998 


0.286 






4000 


0.9985 


0.578 



Table 2: Shapiro- Wilk test under Condition |4] 

Let us now focus on the more general Condition [3] Figure 2 represents the 
empirical distribution of (L/4'\/2 log L) (S^ — ao) in case ao = 3, k = 1 and the 
corresponding narrow-band estimates, whose results are summarized in Table 
3. The improvement in the bias factor with the latter procedure is immediately 
evident. 
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Full Band Estimates 




2 



Narrow Band Estimates 



^210 1 23 



Figure 2: Comparison among biased and narrow estimates (k — 1, L — 2000, 
ao = 3), under Condition |3] 





Shapiro- Wilk Test 


L 


Li 


mean 


var 


W 


p — value 


2000 


1550 


0.072 


0.959 


0.9985 


0.950 


1700 


0.018 


0.951 


0.997 


0.495 


1850 


-0.016 1.004 


0.9977 


0.739 


3000 


2400 


0.092 1.130 


0.9949 


0.920 


2600 


0.072 


0.928 


0.9951 


0.745 


2800 


-0.02 


1.06 


0.9965 


0.340 


4000 


3250 


0.006 


0.985 


0.9968 


0.443 


3500 


0.004 1.097 


0.998 


0.834 


3750 


0.0007 


1.073 


0.9982 


0.874 



Table 4: Normalized Narrow bands data, under Condition |3j k = 1, 

Go = 2. 



ao 



Considering the correction term from Remark |17[ the sample bias is 
consistent with the asymptotic value to three decimal digits. 

Finally, in Figure 3, we report the results obtained on a set of simulations 
under Condition [2] where we have: 



G{1) = Go 1 -f y - ^ 



with Go = 2, ao = 4, L = 4000, Li = 3750. 

We obtain a mean value E (S^ — ao) = 0.040 and a normalized variance 
of 0.9918. Shapiro- Wilk Gaussianity test gives as result W = 0.9981 with a 
p — value = 0.8669. 
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Narrow band estimates 




Figure 3: Distribution of normalized (a^ — Oq) under Condition [2j 

8 Conclusions 

We view this paper as a first contribution in area wliicli deserves much further re- 
search, i.e. the investigation of asymptotic properties for parametric estimators 
on a single realization of an isotropic random field on the sphere. As mentioned 
earlier, an enormous amount of applied papers have focussed on this issue, es- 
pecially in a Cosmological framework, but no rigorous results seem currently 
available. Our results suggest that consistency and asymptotic Gaussianity are 
feasible for spectral index estimators, the rate of convergence being L/logL; 
these estimates are centred on zero in "parametric" circumstances, i.e. where 
the correct model being provided for C/ up to a factor o(j); when the latter 
assumption fails, alternatively, narrow-band estimates can be entertained; these 
estimates ensure convergence to a zero-mean Gaussian distribution, convergence 
rate becoming L/log^ L. 

Many questions are left open by these preliminary results. The first we men- 
tion is the characterization of a whole class of parameters for which asymptotic 
Gaussianity and consistency may continue to hold. More challenging is the 
possibility to relax the Gaussian assumption and consider more general, finite- 
variance isotropic Gaussian fields. In this respect, results in [53] suggest that 
the Gaussianity assumption may indeed play a crucial role, as high-frequency 
consistency and Gaussianity seem very tightly related, for instance when con- 
sidering the asymptotic behavior of the angular power spectrum. It seems also 
important to explore the connection between the spherical estimates we have 
been considering and fixed-domain asymptotic results for Matern-type covari- 
ances which have been discussed on by [33], [I], [2S] and others. 

More interestingly, we aim at relaxing some of the assumptions introduced in 
this paper to make this work more directly applicable on existing datasets. The 
harmonic estimates we have been focussing on hence require the observation of 
the random field on the full sphere. This condition often fails in practice: for 
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instance, in a Cosmological framework large regions of the sky are not masked 
by Foreground sources such as the Milky Way. In ongoing research (see [TU]). 
we are hence considering a Whittle-type estimator based on spherical wavelets 
(needlets, see [55], |27j, [3]), rather than standard Fourier analysis. We believe 
that the material in the present paper can provide also the right starting point 
for these kind of further developments. 



A Appendix 

We start with the proof of Proposition [M] For the sake of brevity, we prove 



only (18), omitting the proof of (16), considering it a particular case in which 



Li — 1. The following Lemma is pivotal for the proof. 
Lemma 21 Let Li < L , then we have 



1 



1 



2(l+i)_^2(l+i) 



(30) 



' log xdx — 



(l+i)_^2(l+i) ^2(l+f) 



L'y'+i)\ogL-Ll^^^^-hogLi 



2 1 



+ 



(31) 



2x^+'' log^ xdx 



L"l^+5JiogL-ii^'^*hogLi 



2(1 + !)^ 



i2(i+f)log2i-Li('+*hog2Li 



+ 

Proof. First of all, observe that for u > 



{21 + 1)1' (log 0" = I 2x^+' log" xdx + OL (1) 



Simple calculations lead to: 



L a;2(l + f) 

2x^+'dx = -J ^ 

(1 + i) 



X — h\ 



L 



!(l+i)_^2(l+i) 



By the same procedure we obtain (31 ) and (32) 



Proof. (Proposition 14 1 We start by observing that 



5](2; + i);nog2A ( ^(2; + i)H 

\i=Li / \;=Li / 



(32) 
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2 



(1 + !)' V2(i + i)' (1 + f) 

^2(l+i)_^2(l+i)- 



+ ^ } ^^1^'^*^ log (1 - g (L)) X 

(1 + i) 

^ 21ogL-log'(l-ff(L)) 1 +Oi(l) ; 



( ^(2Z + l)PlogZ 



l=Li 



2(l+j)_^2(l+i) 



+ V4(i + i)' (1 + i) 



+ +0.(1), 



so we obtain 



r2 



^ ^ . X^(^+i)L^(^+i)V(l-g(L)) 

"-^^^^^^ = 1(171? oTi)^ ""^'^ 



(^+^)((l-(l-g(^)f(^^^))) 

4 (1 + i)' 

L^(^+i)(l-g(L)f(^+t)log^(l-g(L)) 



Observe that 
lo: 



'g\l -g{L)) = (^-g{L)-^9^iL)-^9HL)+0{9'iL))y 
= g^L)+g'{L)+(^^y^{L) + o{g^{L)) , 
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while 

(l-g(L)f(^+*) 1 ^^^^^ , / / , s 



(1 + f) (1 + ^) ' V V 2 



-2g{L)+{2{l + -\-l\g' (L) 



Thus 



4(i + i)' 



l^[l + is + l)g{L) + 



[1 ■ 2. 



+ ^(.'^ + l)(^« + l).9'(i)]+o(iVW) , 
while simple calculations lead to 

L<i+*)(l-.g(L)f(^+5)V(l-.gW) 



(1 + f)' 



^1 ' 2J 

23 1 

Y ^ 24* ~ 12 



By using (18), we have 



7-4(1+1 )„4 /rx 

Zlml) (-) = '-ri^K (s) + o {lS' (L)) , 

(1 + f) 

as we claimed. ■ 

Now we present the result about the analysis of the fourth-order cumulants. 
Lemma 22 Let Ai and Bi be defined as in\2^ and\27^ As L ^ oo, 

^cumlY.(Ai,+Bi,),Y,iAh+Bi,),Y,iAi,+Bi,),Y, i^u + i?u ) U Ol ( ) 

\ ll I2 ^3 ^4 J 

Proof. It is readily checked that 

I Ci Ci Ci Ci I ^,,-3\ 
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cum 



J GLjap) GLjap) GLjap) Gl(qo) [ 
1 Go ' Go ' Go ' Go J 



L8Z^^ ■ ^ [Q'Q'a'Gi 

The proof can be divided into 5 cases: 
1. 



2. 



v. li I2 I3 I4 ) 

= o(^X:P' + i)'i°e'i) = o(!2^); 

V ^2 ^3 ^4 ^ 

1 /y^^o; , im \ \Gl{plq) Gl{oiq) Gl{plq) Gl{oiq) 



3. 



v. l\ I2 I3 I4 ) 

h - ^) -'I i:(-^H-:, { |, efsi, efsi 



^ 1^(2/1 + 1) log Zi I 5^(2/2 + 1) log Z2 



V t3 tS t5 
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V ^1 I2 I3 I4 ) 

^(2Z+1)^ log^ Icum 1 1, |, ^(2/3 + 1) E(2«4 + 1) log | 



G G Q 



= is {E(2^ + 1) log j E(2« + 1)' log' ^-^^ IC{G{G{C^ 

= $ {E(2^ + i)i°g j E(2' + 1) i°g' « = o(^) ; 

V. l\ I2 I3 Ia ^ 

= ^ E(2^ H- W /-m III I ^(2^, + 1) log I 
^ |^(2«, + l)log«,| ^(2/+l)3log3/ct.m ||, |, |,^(2Z, + 

V /l ) I V ^2 ^ 

= {E(2«i + 1) log h I E(2^ + 1)' V /-m 1 1. 1, |, | } 
= |E(2'i + 1) log h I E(2^ + 1) log' ^ = O(^) . 
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